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Abstract: A detailed study of various two-particle correlation functions involving 
photons and neutral pions is presented in proton-proton and lead-lead collisions at 
the LHC energy. The aim is to use these correlation functions to quantify the effect of 
the medium (in lead-lead collisions) on the jet decay properties. The calculations are 
carried out at the leading order in QCD but the next-to-leading order corrections are 
also discussed. The competition between different production mechanisms makes the 
connection between the jet energy loss spectrum and the 7—71° correlations somewhat 
indirect while the 7 — 7 correlations have a clearer relation to the jet fragmentation 
properties. 
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1. Introduction 

Electromagnetic probes have long been thought to be useful to detect the formation 
of a quark-gluon plasma in ultrarelativistic heavy ion collisions [1, 2, 3, 4]. Many 
observables involving photons can in principle be used. One of the simplest, from a 
theoretical point of view, is the single photon spectrum as a function of the transverse 
momentum p T : it is expected that secondary collisions in heavy ion scattering will 
produce an excess of direct photons, as compared to proton-proton scattering, in an 
energy domain a few times the plasma temperature. However the flux of background 
photons from hadronic resonances is quite large and this makes the extraction of 
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the direct photon signal a non trivial task experimentally. These observables would 
probe the plasma at the early times of the collision when the medium is the hottest. 



Photons can also be used in a different kinematical regime, namely at large 
transverse momentum (many times the temperature of the plasma). Such photons 
are produced in primary collisions i.e. as in proton-proton collisions and, ideally, 
their production rate is calculable in perturbative QCD in the next-to-leading order 
(NLO) approximation [5]. For a high enough transverse momentum these direct 
photons should be easily extracted from the background since the ratio 7/71° is rapidly 
increasing with p T . One can then study the decay properties of the jet recoiling 
from the photon by considering various photon-hadron or photon-photon correlation 
functions where the hadron or the second photon are fragments of the jet [6, 7, 8]. In 
the simplest case, when the transverse momentum of the recoiling jet exactly balances 
that of the photon, such observables allow to map out the fragmention function of the 
jet traversing the medium [9]. For a sample of data large enough one hopes to study 
the difference in the shape of the fragmentation functions in proton-proton and ion- 
ion collisions. Of course the real situation is more complicated because the large p T 
photon can itself be produced by bremsstrahlung [10] in which case the photon and 
the jet momenta become somewhat uncorrelated as it is also the case when higher 
order corrections are taken into account. Furthermore, in order to have a reasonable 
counting rate for the correlation studies one cannot consider photons with too large 
transverse momenta and the 7T° background may then remain a problem. Thus if 
one studies 7 — n° correlations, the tt° — ir° contribution should also be considered 
in turn. 



In the following we will study various 7 — 7 , 7 — 77° and n° — 7r° correlations both 
in proton-proton and lead-lead collisions at the LHC. Shadowing effects will be con- 
sidered and, following standard practice, we will assume that the effects on the hard 
process of the parton multiple scattering through the medium can be parametrized 
by a modification of the fragmentation functions. The results are obtained in the 
leading-logarithm approximation of QCD since the status of NLO calculations in a 
medium is not yet clear. Our results should therefore be considered only as semi- 
quantitative. The model will be presented in the next section with special emphasis 
on the medium modified fragmentation effects. Then we turn to several observables 
and compare their behavior in proton-proton and lead-lead collisions. A discussion 
of the effects of NLO corrections in proton-proton collisions is given specifically to 
test the stability of the shape of observables. We consider this to be indicative of 
the stability of correlations under higher order corrections in heavy ion collisions. 
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2. The model 



2.1 The leading order cross section 

At the leading order in QCD the basic two-particle correlation cross-section, from 
which we will construct various observables, can be written [6] 

da AB ^ CD 1 ^ D c/c (z 3 ,M F )D D/d (z 4 ,M F ) 

dp T3 dy 3 dz 3 dp T4 d yi dz 4 8ns 2 ^ a > b > c ' d z 3 z 4 1 T3 tJ 

F a/A (x u M) F b/B {x 2) M) 

x, x~ 2 {2A} 

where p Ti and y iy i — 3, 4, are respectively the transverse momenta and rapidities of 
the final state particles. The momentum k Ti , % = 3, 4, is the transverse momentum of 
the parton c (respectively, d) which emits particle C (respectively, D) of momentum 
p Ti , carrying a fraction Zi = p Ti /k Ti of the parton momentum. The fragmentation 
functions D c / C and D D / d depend on the collinear factorization scale M F . For the 
production of hadrons we will use the leading order functions of [11]. When a pho- 
ton is detected in the final state it can be either produced directly, in which case 
the fragmentation function D 1 / c (z,M F ) reduces to a Dirac function 5(1 — z) or it 
can be produced via bremsstrahlung of a final state quark or gluon (see Figure 1 
for illustration). In the latter case we use the BFG parametrization 1 of [12]. The 
structure functions F a /A and F h / B of the projectile and target, A and B, depend on 
the factorization scale M and they are normalized to one nucleon. Our standard in 
the following study is the parametrization of CTEQ6L [13]. 

The quantity \M\ 2 ah ^ cd is the matrix element squared, averaged over spin and 
color, of the partonic sub-process ab — > cd. It depends implicitely on the renormali- 
sation scale \i through the strong coupling a s (fi). In the following, unless otherwise 
specified, all scales, //, M, and M F , are set equal to (p T + p T )/2. Higher order cor- 
rections to Eq. (2.1) have been calculated and extensive phenomenological studies 
have been made in the case of proton-proton collisions [6, 7, 8, 14, 15, 16]. They are 
briefly discussed in Section 6. 



2.2 Initial state nuclear effects 

The structure function of a nucleon in a nucleus is modified by shadowing and anti- 
shadowing effects. These effects are hard to calculate theoretically and in this study 
we will simply use a parametrization of the parton distribution functions measured in 
deep-inelastic scattering experiments of leptons off nuclei [17] and Drell-Yan produc- 
tion in proton- nucleus reactions [18]. We follow the approach of Eskola et al. [19] who 

1 These parametrizations are given at the NLO of QCD. Nevertheless we use them for our leading- 
logarithmic studies for lack of recent leading- logarithmic parametrizations. 
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Figure 1: Two processes which contribute to 7 — it production at leading order: the 
photon may be produced directly (one-fragmentation, labeled If in the text, left) or by 
parton fragmentation (two- fragmentation, labeled 2f, right). 

tabulate a function S a / A {x,M) which relates the parton distributions in a nucleon 
N to those in a nucleus A via 

F a/A (x, M) = S a/A (x, M) F a/N (x, M) (2.2) 

Unlike older parametrizations, it is to be noted that the parametrization of [19] treats 
quarks and gluons separately and, furthermore, the factorization scale dependence 
is taken into account. In this parametrization, shadowing effects reduce the quark 
(by up to 10%) and the gluon (by up to 30%) distributions in the nucleon below 
x < 0.03, while ant i- shadowing enhances the distributions by up to 10% between 2 
x = 0.03 and x = 0.4. Since observables usually involve an integration over a rather 
large range in x, such shadowing/antishadowing effects will not affect the predictions 
very much. A summary of recent shadowing studies at LHC can be found in [20]. 

Plugging Eq. (2.2) in Eq. (2.1), with proper account of isospin effects, one obtains 
the cross section normalized per one nucleon in a nucleus. To obtain the counting 
rate in an actual heavy ion experiment one needs to account for the number of 
nucleon-nucleon scatterings occuring in a nucleus-nucleus collision. This is done in a 
standard way using Glauber theory. A hard cross section for an AA collision, with a 
given centrality class C (equivalently impact parameter range), is obtained from the 
corresponding hard cross section via the "binary scaling" relation 

„hard\ /at \| AA Jiard in q\ 

a AA \C = {^coll)\c &NN ( 2 - 6 ) 

&NN 

where (N co u)\ c is the number of collisions at the chosen centrality, o A eo A is the geometric 
cross section obtained via the Glauber multiple scattering model and <jnn is the 

2 The precise x values depend on the factorization scale. 
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nucleon-nucleon cross section. Details, as well as the numerical values of the various 
terms, are given in [4]. In particular, for collisions with a centrality less than 5% the 
estimate is {N co u)\ c = 1876 and cr^ = 7745 mb for lead-lead collisions at 5.5 TeV 
with <r NN = 72 mb. In the following we implicitely consider only central collisions, 
with C < 5%, for lead-lead collisions. When quoting numbers of events for a given 
observable we assume the standard luminosity for lead-lead collisions in ALICE, 
£ = 5.10 26 cm -2 sec -1 , and make the hypothesis that LHC is running 30 days per 
year in the heavy ion mode. 

2.3 Medium-modified fragmentation functions 

Much progress has been done over the last decade to better understand the gluon 
radiation by hard partons travelling through dense QCD media [21]. More recently, 
the important connection between the medium- induced gluon spectra dl/du and 
the probability distribution V in the energy loss has been made explicit [22] and 
computed numerically soon after [23, 24] . However, it remains unclear how to relate 
the parton energy loss mechanism to observable quantities. 

For sufficiently large k T parton production, nevertheless, a clear separation is 
achieved between the hard production process, with a time scale O (& T _1 ) , the effects 
of the medium, O (t me d), and the fragmentation mechanism, O (k T /A 2 ), 

^- « tmcd « ^ (2.4) 

Provided the hierarchy (2.4) is justified, it is sensible to model the energy loss 
effects at the level of fragmentation functions. In the present study, we shall follow 
the model suggested in Ref. [9] in which the energy e lost by the hard parton leads 
to a rescaling of the momentum fraction Zd 

z& — _> z * — ^ Td — Zd (2 5) 

k Td k Td — e 1 — e/k Td 

in presence of a QCD medium. Consequently, the medium-modified fragmentation 
functions D™ d (zd, M F , k Td ) may simply be expressed as a function of the standard 
(vacuum) fragmentation functions D D /d{zd, Q 2 ), 

z d D™ / d d (z d ,M F ,k Td )= / deV d (e,k Td ) z* d D D/d (z* d , M p ). (2.6) 

Jo 

Here, Pd(e, k Td ) denotes the probability for the parton with energy k Td to lose the 
energy e [22], which has been given a simple analytic parametrization in [23] which 
we shall use in the present calculations. The calculation of Ref. [23] is based on the 
medium-induced gluon spectrum determined by Baier, Dokshitzer, Mueller, Peigne 
and Schiff (BDMPS) including 0{l/k Td ) corrections [25, 26]. Such a model was 
shown to describe successfully hadron production in semi-inclusive DIS reactions off 
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nuclear targets [27]. The BDMPS framework should be particularly suited when the 
number of scatterings incurred by the hard parton in the QCD medium (opacity) is 
large. While thick and dense media are indeed expected to be produced in nuclear 
collisions at LHC energy, we note however that such a calculation may not properly 
describe the energy loss process for partons produced close to the surface. Let us 
mention that the probability distribution at finite opacity was also determined in the 
soft limit (k Ti — > oo) in Ref. [24] which lead to recent phenomenological applications 
at RHIC and LHC [28, 29]. 

The BDMPS energy loss distribution is characterized by the energy scale [25] 

co c = \qL 2 (2.7) 

where the so-called gluon transport coefficient q reflects the medium gluon den- 
sity [30] and L the length of matter covered by the hard parton in the medium. Note 
that q in (2.7) has to be seen as a time averaged quantity (q) to take properly into 
account the longitudinal expansion of the produced medium [31]. 

For the calculations to come, a qualitative estimate of u c for the dense medium 
produced in lead-lead collisions at LHC energy is needed. The transport coefficient q 
is directly related to the gluon density whose increase from RHIC to LHC is of order 
6 — 7 in hydrodynamical calculations [28, 32], while a smaller increase 2 — 4 in the 
hadron multiplicity at mid-rapidity (also linked to q) is predicted by several models 
(a review can be found in [33]). Using the estimate based on the pion p T spectra 
measurements at RHIC, ^c| RHIC — 10 — 20 GeV, within the same framework [23, 34] 
we shall take throughout this study the rather conservative choice 3 , u c = 50 GeV. 

Let us remind the reader that the goal here is not to provide quantitative predic- 
tions but rather to show typical trends one could expect in 7—7 and 7— 7r° correlations 
in heavy- ion collisions at the LHC. Therefore, our conclusions should not depend 
much on the precise value we assume for the energy loss parameter u c . 

Since fragmentation functions fall steeply with z, even a small shift Az^ = z* d — 
za ~ Zdt/k Ti in Eq. (2.5) may substantially affect the fragmentation process due to 
parton energy loss. This can be seen for instance in Figure 2 where the fragmentation 
functions into a photon and into a neutral pion, using respectively the BFG [12] and 
KKP LO [11] parametrizations, are computed for k T = 50 GeV up quark and gluon 
traversing the medium (u c = 25, 50 GeV) or not {uo c = GeV). 

First, Figure 2 indicates that medium effects prove stronger for gluon than for 
quark fragmentation. The origin is actually twofold. First, hard gluons lose more 
energy than quarks do from their larger color charge (C g — 3, C q — 4/3). Moreover, 
the quenching of medium-modified fragmentation functions Eq. (2.6) increases with 

3 This choice is motivated by the fact that the gluon distribution at very small x could evolve 
more slowly than seen so far at HERA due to possible gluon saturation. Moreover, we want our 
predictions to be seen as lower estimates as far as medium effects are concerned. 
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Figure 2: Medium-modified fragmentation functions D^J^Zd, M F , k Td ) for various energy 
loss scales, oj c = (vacuum), 25 and 50 GeV. The parton energy is k T = 50 GeV and the 
fragmentation scale is set to M F = p T /2. 



the slope of vacuum fragmentation functions, much steeper in the gluon channel. 
Finally, we observe that the effects of parton energy loss become more pronounced 
as z gets larger, due to the restricted available phase space in Eq. (2.6). 

The medium-modified fragmentation functions depend now explicitely on the 
parton energy, k T . To show the sensitivity of the medium effects on the parton energy, 
the ratio of medium (using uj c = 50 GeV) over vacuum fragmentation functions is 
determined for k T = 25, 50 and 100 GeV quarks and gluons. As can be seen in 
Figure 3, medium effects will be magnified as the parton energy is getting smaller. 
When k T becomes too small as compared to uj c , however, the picture of a hard 
parton penetrating the soft medium is no longer correct and the applicability of the 
energy loss framework becomes doubtful. It was shown for instance in [23] that the 
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k T = 25GeV — k T = 50GeV 



k T = 100 GeV 
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Figure 3: Ratio of medium-modified (u> c = 50 GeV) over vacuum (o; c = GeV) fragmen- 
tation functions for various parton energy, k T = 25, 50 and 100 GeV. The fragmentation 
scale is set to M F = p T /2. 

eikonal approximation is explicitely broken down for parton energy smaller than half 
the scale u c . This has to be seen as the lower limit for the most energetic photon 
transverse momentum. On the other hand, in the high energy limit k T 3> uo c and 
thus z* ~ z, the medium effects vanish and the ratio approaches one. 



3. The correlations 
3.1 Observables 

From the kinematical variables available in Eq. (2.1) we can construct the following 
observables: 
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- the invariant mass of the particle pair, 



m u = 2 (PtzPt^HVS ~ 2/4) - P T z ■ Pta) 

= 2z 3 z 4 kl (ch(y 3 -y4) + l) (3.1) 



the transverse momentum of the pair, 



\PtZ PtA I 

k T \z 3 — z A \ (3.2) 



- the relative transverse momentum of the particles (also called momentum bal- 
ance [6, 35]) 

Ptz ■ Pt4 
z 34 — n 

= Z f (3.3) 

where k T is the common value of the transverse momentum of the final state partons. 
For completeness we quote here the expressions of the Xi values of the initial partons, 

Xl = -L ( + ^e yi 

X2 = JL ( Eil e -V3 + Pl± e -vA (3.4) 
V s V z 3 z A J 

In the case of 7 — n° correlations, one has z 3 = 1 when the photon is produced 
directly. Fixing furthermore the rapidity of the photon and the pion in a narrow 
range around 0, for example, we are left with two independent kinematical variables 
z A and k T and the expressions of the observables defined above considerably simplify. 
One has: 

ml 4 = Ak 2 T z A 

Qrp — — krp J l Z A J 

Z34 = z A (3-5) 

which show a straightforward relation between the fragmentation variable z 4 and the 
observables. 

When studying the observables we integrate over k T above a given value. Since 
the cross section is rapidly falling when k T is increasing the effective transverse 
momentum will remain close to its minimum value leaving z 4 as the only effective 
variable. From the behavior of the above observables one should get constraints on 
the behavior of the fragmentation function if we assume that the structure functions 
are precisely known. 
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When the photon is produced via bremsstrahlung or in the case of 7r°— 7r° correlations 
the above simple situation is somewhat smeared because Z3 is now a relevant kine- 
matical variable. However, when studying asymmetric configurations with a large 
p T3 particle on one side and a small p T4 particle on the other side, trigger bias ef- 
fects will force a large value of z 3 and some correlation is still expected between the 
observables of Eqs. (3.1) to (3.3) and the z 4 dependent fragmentation function. 

On the other hand, considering 7 — 7 correlations, when both photons are pro- 
duced directly, one has an over constrained system with Z3 = z 4 = 1 and the dis- 
tributions in Z34 and in q T reduce to Dirac 5 functions while the invariant mass 
7?7| 4 = 2/^(ch(7/3 — 7/4) + 1) is regular. It is obvious that higher order corrections will 
smear the 5-function singularity. More precisely, Z34 = 1, or equivalently q T = 0, is 
an infrared sensitive point and an accurate prediction of the behavior of these ob- 
servables near this point will require the resummation of large \n 2 (q^/s) and ln(^/s) 
terms. 

3.2 Kinematical cuts 

We study a basic perturbative QCD (pQCD) mechanism modified by the presence of 
a dense environment. It is necessary to insure that the particles we observe are decay 
fragments of jets and are not produced by secondary collisions. Recent studies [4] 
in the framework of perturbative QCD for primary collisions and a hydrodynamic 
model to describe secondary collisions have shown that, at LHC, particles produced 
above p T = 5 GeV are of pQCD origin. We therefore impose a minimum transverse 
momentum of 5 GeV on the particles from which we construct the various correlation 
observables. 

To study a large domain in the fragmentation variable z it is necessary to consider 
asymmetric configurations. Another constraint is to be able to distinguish photons 
from pions which requires, for ALICE for example, p T ^> 25 GeV. On the other 
hand, to have a reasonable counting rate, one should not go to too high values of 
p^ . Besides, if jets are too energetic, energy loss effects will be small and difficult 
to observe. 

In the subsequent studies, apart from various p T distributions, we will look at the 
three distributions in Eqs. (3.1) to (3.3) with p T > 25 GeV and p T ^ > 5 GeV to satisfy 
the above criteria. We will also consider higher cuts in p T to probe the sensitivity of 
the energy loss mechanism on the jet energy. When displaying the distributions we 
always assume photons and pions are produced in an interval of 8y — ± 0.5 unit of 
rapidity around y — 0. 

4. Phenomenology of 7 — 7r° correlations 
4.1 Dynamical components 

Before studying the shape of correlation functions in lead-lead collisions it is worth- 
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while considering into some details the case of proton-proton scattering to better 
understand the dynamics of the reaction. As above mentioned (see Figure 1), the 
photon can be produced directly and only the recoiling jet fragments into a pion 
(labeled If, open squares in the following figures), or both the photon and the pion 
are produced by fragmentation of partons (labeled 2f, full squares). 
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Figure 4: The If and 2f components for various kinematical configurations in 7 — 
7T° production for proton-proton scattering at yfs = 5.5 TeV. Top: The p T ^and the 
p T7r distributions for the cuts p T ^ > 25 GeV and p T ^ > 5 GeV. Bottom: Same as above 
for p T > 50 GeV and p T > 5 GeV. 



The relative weight of each mechanism depends crucially on the transverse mo- 
mentum of the produced particles. This is illustrated in Figure 4. On the right 
panels one displays the production rate as a function of p T ^ : typically, the 2f process 
dominates when p T ^ < 50 GeV but decreases much faster than the If component as 
p T increases and becomes negligible for photon transverse momenta above 100 GeV. 
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Therefore, changing the p^cut will affect dramatically the If and 2f relative contri- 
butions to 7 — 7T° production. 

As shown in Figure 4 (left), the pion p spectrum is dominated by the 2f process 
in a large p T ^ range when the photon energy is not too large (p T > 25 GeV) while 
applying a higher cut (p T > 50 GeV) in p T ^ results into a comparable magnitude 
for both mechanisms on a wide p T ^ domain. Unlike the p T distribution the relative 
weight of the 2f process is increasing with the transverse momentum of the pion in 
the kinematical domain shown. Indeed, the If component is disfavored at large p Tir as 
it requires more energetic photons {p Ti > P T7T ), while the 2f process allows the photon 
to keep a small transverse momentum, slightly above the 25 GeV or 50 GeV cut. 

Experimentally, these two contributions may be disentangled by means of calorime- 
try techniques using appropriate isolation criteria. However, the large multiplicity 
reached in high energy heavy-ion collisions prevent one from using such techniques. 
It should also be reminded that when performing a full NLO study, the distinction 
between the leading-order fragmentation and the next-to-leading order direct com- 
ponent is arbitrary and depends on the fragmentation scale Mp, only the sum of 
these components being meaningful and ideally scale independent [7]. 

Medium effects may change considerably whether one or the other process dom- 
inates since the typical parton energy, k T = p T /z 3 , is quite different for If (z 3 = 1) 
and 2f (z 3 < 1). Naively, the effects of parton energy loss should be stronger when 
both the pion and the photon come from the hard parton fragmentation. However, 
this is not necessarily true since the parton energy is actually much greater in the 2f 
channel, for which medium effects prove weaker (see end of Section 2.3). We shall 
come back to these observations when discussing the proton-proton and lead-lead 
spectra in the next section. 

4.2 Distributions 

In Figure 5 we discuss four distributions, respectively in the pion transverse momen- 
tum p T ^ , the photon transverse momentum p T ^ , the 7 — ir° invariant mass m vr7 and 
the transverse momentum of the pair q T . We impose the following cuts: p Tjr > 5 GeV 
and p T > 25 GeV. In each case three curves are displayed: proton-proton scatter- 
ing (open dots), lead- lead scattering with shadowing but without energy loss (full 
squares) and lead-lead scattering with shadowing and energy loss using ui c = 50 GeV 
(open squares). 

It is clear that shadowing and isospin effects do not modify the distributions very 
much: a small antishadowing effect can be observed at large transverse momenta or 
at high invariant mass due to the fact that the kinematics then becomes sensitive to 
larger x partons in the nuclei. Energy loss effects are quite visible particularly at the 
low p T values of the pion or the photon. On the other hand, to produce a pion at 
high transverse momentum requires a parton with large k T for which the energy loss 
is expected to be smaller. We observe, accordingly, that the spectrum in lead-lead 
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Figure 5: The four distributions in 7 — tt° production defined in the text for proton-proton 
(open dots) and lead-lead scattering (no energy loss: black squares; with energy loss: open 
squares) at = 5.5 TeV. Both the photon and the pion are produced at rapidity [-0.5, 
0.5] and the following cuts are imposed: p T > 25 GeV and p T ^ > 5 GeV. 



collision tends to approach the proton-proton spectrum as p T increases. The medium 
effects are also particularly visible on the spectrum as a function p T< : as long as the 
photon is produced directly (If), the p T ^ spectrum reflects the energy of the parton, 
k T = p T ^ , which eventually fragments into the pion. Again, the quenching will be 
maximal for small p T ^ (small k T ) while at asymptotic energies, parton energy loss 
will have no observable consequence. Similar behavior is observed in the invariant 
mass distribution: small masses correspond to low k T partons and therefore lead to 
a stronger suppression. One may notice, in passing, the rather large counting rates: 
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with the numbers given in Section 2.2 we estimate that about 1300 7 — 7r° pairs with 
a 100 GeV invariant mass will be produced, per year, in ALICE. 

Perhaps, more interesting is the q T spectrum which exhibits a maximum when 
the pion and the photon transverse momenta lie just above the imposed kinematic 
threshold, which is located at the difference between the p T and the p T ^cut, 20 GeV. 
Above 20 GeV, the distribution is reminiscent of the p T and the m nj distribution 
and, in particular, the larger the q T the weaker the energy loss effect. Similarly, the 
energy loss effects will tend to be smaller at very small q T <C 20 GeV as the pion 
transverse momentum and thus the 2f contribution - less affected by the medium - 
increases with decreasing q T . Therefore, we expect the medium effects to be maximal 
for q T roughly around the difference of the transverse momentum cuts. 

All these features are best seen when normalizing, to the proton-proton distribu- 
tions, the lead-lead distributions with energy loss (open squares) or without energy 
loss (full squares) in Figure 6. In all cases it appears that the observables are affected 
by antishadowing, and not shadowing, but this remains a small effect, less than 10% 
in general. Energy loss effects, on the other hand, modify the distributions much 
more drastically. The p T ^ spectrum is suppressed by about 50% below p TiT = 25 GeV 
but no suppression occurs above 50 GeV. As for the p T ^ distribution, the suppres- 
sion is maximum at low transverse but is monotonously reduced as the momentum 
increases. In the q T spectrum the change of slope discussed above it is particularly 
noticeable. 

The same spectra are computed in Figure 7 assuming a larger cut for the pho- 
ton transverse momentum, p T > 50 GeV. Although the If contribution becomes 
relatively more important, one observes similar features as before. Again, the nor- 
malized q T spectrum shows a clear minimum in Figure 8 around q T ~ 30 GeV, under 
which the 2f contribution starts to dominate. The quenching of these spectra proves 
less pronounced - the ratio decreases down to 0.6 in the q T spectra to be compared 
to 0.35 before - since the initial parton energy is twice as large. Finally, let us note 
that the counting rate drops by a factor 5 to 10 when increasing the photon cut from 
25 GeV to 50 GeV. 

We turn now to a detailed discussion of the distribution in the momentum bal- 
ance z 34 for both cuts p T ^ > 25 GeV and p T ^ > 50 GeV, keeping p T ^ > 5 GeV as before 
(Figure 9, top). The maximum of these distributions is reached for z 34 = 0.2 and 
0.1, i.e. the ratio of the pion over the photon transverse momentum cuts. Smaller 
(larger) z 34 values are obtained by increasing the photon (pion) transverse momen- 
tum. Looking at the medium effects (Figure 9, bottom), one observes rather struc- 
tureless features: as compared to the proton-proton case the spectrum is reduced 
to 40% (respectively 60%) over most of the z 34 range when the cut p T > 25 GeV 
(respectively p Ti > 50 GeV) is imposed. Below z 34 < 0.2, the suppression is not as 
strong since the photon energy, hence k T , is getting larger. 

We stressed in Section 3 that the momentum balance z 34 is closely related to the 
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Figure 6: 

tering. 



Same as Figure 5 but the distributions are normalized to proton-proton scat- 



fragmentation variable z when the photon is produced directly. Therefore, it may 
look surprising at first glance not to observe the suppression becoming more impor- 
tant when z 34 gets close to one, as the ratio of medium over vacuum fragmentation 
functions may suggest in Figure 3. To understand the origin of the flat behavior of 
this ratio, the individual If and 2f contributions to the z 3i distribution are repre- 
sented in Figure 10 (top). When the p T ^cnt is set to 25 GeV, most of the events (for 
z 34 > 0.15) actually come from the double fragmentation process and the relative pro- 
portion of this 2f contribution increases with z 34 . For this 2f component, however, the 
connection between the momentum balance and the fragmentation variable is lost. 
In particular, the pion transverse momentum gets larger as z 34 increases, thereby 
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Figure 7 : Same as Figure 5 with the following cuts: p T ^ > 50 GeV and p T ^ > 5 GeV. 

reducing the medium effects. To be more explicit, we show in the lower panel of 
Figure 10 the medium effects on the individual If (open squares) and 2f components 
(full squares) separately. As expected, the suppression in the 2f channel becomes less 
important when z 34 increases unlike the If channel, whose suppression is reminiscent 
to the z dependence of the medium over vacuum fragmentation functions, with a 
vanishing ratio at z 3i ~ 1. Summing the two contributions, the resulting suppression 
(open circles) is an interplay between the If and 2f behavior. As far as counting 
rates are concerned they are rather large: for ALICE one expects about 5 10 5 pairs 
at z 34 = 0.5 when p T ^> 25 GeV and 5 10 4 pairs when p T ^> 50 GeV. 

In order to isolate the If channel - whose medium effect is remarkable - it 
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Figure 8: Same as Figure 6 with the following cuts: p T > 50 GeV and p T > 5 GeV. 



would be necessary to increase the photon transverse momentum, making the 2f 
process highly unlikely. It can be seen from the right panel of Figure 10 that going 
from a 25 GeV to a 50 GeV p T ^cut indeed increases significantly the If component. 
Nevertheless, the 2f contribution remains too large at large z 34 to observe a huge 
medium suppression in this kinematical region. It may then be necessary to trigger on 
even more energetic photons, the drawback of too stringent cuts being the smallness 
of the corresponding cross sections. 

As we shall see in the next Section, diphoton production mostly comes from the If 
contribution process at the LHC. This observable may therefore be more interesting 
than 7 — 7T° correlations, at least regarding the momentum imbalance distributions. 
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Figure 9: The z 34 distribution in 7 — 7r° production for proton-proton (open dots) and 
lead-lead scattering (no energy loss: black squares; with energy loss: open squares) at 
y/s = 5.5 TeV. Both the photon and the pion are produced at rapidity [-0.5, 0.5] and the 
following cuts are imposed: p T ^> 25 GeV and p Tir > 5 GeV (left) and p T ^> 50 GeV and 
p Tn > 5 GeV (right). Bottom: The same distributions normalized to the proton-proton 
case. 



5. Phenomenology of 7 — 7 correlations 
5.1 Dynamical components 

This section is devoted to the study of 7 — 7 (or diphoton) correlations in the same 
kinematic regime as before. As compared to the previous cross sections they are, 
in principle, a factor 0(a/a s ) smaller, but the counting rates at the LHC should 
nevertheless remain sufficient for our studies. On top of the If and 2f components, the 
new feature is that both photons can be produced directly (direct process, Figure 11) 
in which case they are not affected by the medium. 
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Figure 10: The If and 2f components in 7 — 7r° production for various kinematical con- 
figurations for lead-lead scattering at y/s = 5.5 TeV. Top: The z 34 distribution for p T ^ > 
25 GeV (left) and p T ^ > 50 GeV (right). Bottom: The same distributions normalized to the 
proton-proton case. 



The relative weight of all three components, for proton-proton scattering, are 
shown in Figure 12 imposing a minimum transverse momentum of 25 GeV on one 
photon and 5 GeV on the other. We first consider (top left) the spectrum in p T ^. It 
shows a discontinuity at 25 GeV since below this value only one photon is measured 
and, furthermore, only the If and 2f processes contribute. One observes the domi- 
nance of the If component for the whole transverse momentum range, even when the 
direct process contributes, i.e above the larger p T ^cut. Concerning the invariant mass 
spectrum, the 2f component is at least three times smaller than the If contribution, 
while the direct piece causes a small hump to the total cross section at twice the 
p T threshold value. Note the "singular" contributions of the direct component at 
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Figure 11: On top of the one-fragmentation and two-fragmentation process, two photons 
can be produced directly in 7 — 7 production at leading order. This process is not affected 
by the medium and yields singular contributions at q T = GeV and z 34 = 1. 

q T = and z 34 — 1, a feature of the leading logarithmic approximation. The ordinate 
of the corresponding points depend, of course, on the size of the bin (equivalently the 
resolution of the detector). The shape of the distributions near this infrared singular 
point is expected to be modified by the higher order corrections. 

5.2 Distributions 

We now compare in Figure 13 the p T and q T spectra in proton-proton scatter- 
ing (open dots), lead- lead scattering with shadowing but without energy loss (full 
squares) and lead- lead scattering with shadowing and energy loss (open squares). 
The ratios of the lead-lead over proton-proton spectra are also shown. Similarly to 
the 7 — 7T° case, the effects of nuclear shadowing turn out to be negligible. On the 
contrary, interesting features due to the energy loss mechanism are observed. The 
strongest suppression of the p T ^ spectra is reached for transverse momenta of order of 
the upper cut, p T ^ 20 GeV. This can be understood as follows. As p T ^ approaches 
the upper cut from "below" , events with larger z are selected, p T ~ p T , where 
energy loss effects are most pronounced (cf. Figure 2). Above that cut, the pro- 
portion of directly produced photons (unaffected by the medium) is getting larger 
and the quenching factor is slowly reaching unity as expected at asymptotic ener- 
gies, p Tni ^> u c . Looking at the ratio of q T spectra, the smaller the q T the larger the 
suppression of diphoton events in lead— lead collisions. Indeed, since the production 
is dominated by the If process, with 23 = 1, small q T events imply a large value 
for the other photon (q T = k T \ Z3 — z±\) where medium effects are the strongest. 
Moreover, we no longer observe the same feature as in 7 — n° production - the mild 
increase of the ratio at very small q T (Figure 5, lower right) - since the 2f fragmen- 
tation contribution to diphoton production is much smaller. Finally, the ratio at 
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Figure 12: The three components (direct, If and 2f) contributing to diphoton production 
in proton-proton collisions at = 5.5 TeV. A lower cut of 25 GeV is imposed on the 
transverse momentum of one photon and 5 GeV on the other. The four distributions 
shown are in the transverse momentum of either photon, the diphoton invariant mass, 
the transverse momentum of the pair, the diphoton invariant mass and the momentum 
imbalance. 

q T = GeV is almost close to one due to the singular contribution of the direct 
process, unaffected by the medium. 

The distribution in the diphoton momentum imbalance z 34 is now discussed in 
Figure 14, using the previous kinematical cuts p T ^ > 5 GeV and p T ^ > 25 GeV 
(left) as well as p T ^ > 5 GeV and p T > 50 GeV (right) to keep the parallel 
with Section 4. Similarly to the 7 — 7r° distribution, the distribution is maximal 
around the ratio of the p T ^cuts (z 3i = 0.2 and 0.1 respectively) and decreases rapidly 
with z 34 , a shape reminiscent of the photon fragmentation functions in Figure 2. 
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Figure 13: Top: The photon p T ^and photon pair q T transverse momentum distributions 
in 7 — 7 production for proton-proton (open dots) and lead-lead scattering (no energy 
loss: black squares; with energy loss: open squares) at \fs = 5.5 TeV. Both photons are 
produced at rapidity [-0.5, 0.5] and the following cuts are imposed: p Tll > 25 GeV and 
Pt-,2 ^ ^ GeV. Bottom: The same distributions normalized to the proton-proton case. 

In particular, it is remarkable to notice how the diphoton quenching in Figure 14 
[bottom, open squares) ressembles the ratio of medium- modified over vacuum parton 
to photon fragmentation functions (Figure 3). Unlike the 7 — tt° case, the smaller 2f 
contribution to diphoton production does not spoil too much the large z 34 suppression 
and make the interpretation of the momentum imbalance spectra much easier in 
terms of photon fragmentation functions. 

Notice that one expects a reasonable number of photon-photon events: for ex- 
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ample, for z 34 = 0.5 one has 2 10 4 events for p T _> 25 GeV and 2.5 10 3 for p T ^> 
50 GeV with ALICE luminosity for one month running time. 
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Figure 14: The z 34 distribution in 7 — 7 production for proton-proton (open dots) and 
lead-lead scattering (no energy loss: black squares; with energy loss: open squares) at 
■\fs = 5.5 TeV. Both photons produced at rapidity [-0.5, 0.5] and the following cuts are 
imposed: p Tll > 25 GeV and p T ^ 2 > 5 GeV (left) and p Tll > 50 GeV and p Tl2 > 5 GeV 
(right). Bottom: The same distributions normalized to the proton-proton case. 



6. Qualitative effects of NLO corrections 

As already mentioned, the not yet clarified status of NLO QCD calculations in the 
medium has lead us to consider 7 — 71° and 7 — 7 production to leading order ac- 
curacy. Nevertheless, it is instructive to study and quantify the role of higher order 
corrections in proton-proton collisions. In particular, attention should be paid to the 
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Figure 15: Comparison of LO and NLO calculations for various ir — 7 correlations. 



phase space which gets modified by higher order QCD corrections: consequently this 
can affect our leading order predictions in lead-lead collisions. 

The NLO predictions are obtained using the CTEQ6M structure functions [13], 
the KKP next-to-leading logarithmic fragmentation functions for the 7r° and the 
BFG, set II, fragmentation functions into a photon. As before all scales are set equal 
to (p Tg + p T4 )/2. A cut in the azimuthal angle between the two particles has been 
applied, > min = tt/2. 

The 7 — 7T° correlation functions are plotted in Figure 15 at LO (dots) and NLO 
(squares). The overall effect of higher order corrections proves quite small - say less 
than 30 % - except in some specific regions of phase space. As shown in Figure 16, 
the ratio of NLO over LO spectra is almost constant except at small invariant mass 
or momentum imbalance. The reason comes from the new parton configurations in 
momentum space available at higher order. Indeed, the two particles are no longer 
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Figure 16: Ratio of NLO over LO estimates of the tt — 7 correlations shown in Figure 5. 

constrained to have opposite momenta when going from the two-body to the three- 
body NLO kinematics. This will affect, in particular, variables like the invariant 
mass, the pair transverse momentum or the momentum imbalance which now depend 
explicitely on the relative azimuthal angle <fi between the two particles (see Eqs. (3.1) 
to (3.3)). This can be seen in Figure 16 where the momentum imbalance spectrum 
is clearly enhanced at NLO when Z34 oc cos</> gets very small: this corresponds to 
configurations where the observed particles are recoiling from the third undetected 
jet in the opposite hemispere. Other effects include the shift of the threshold in 
the invariant mass distribution, from m vi = ^J^p Ti P T7i = 10v5 GeV (</> = 7r ) to 
m vi = ^J2p~p^~ = 5-\/T0 GeV (0 = 7r/2), or the shift of the pair momentum 
spectrum to larger q T which results in the ratio smaller than one below 20 GeV (the 
difference of the p T cuts, i.e. the maximum of the distribution) and larger above. 
All these effects depend crucially on the cut m j n in the azimuthal angle and should 
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Figure 17: The NLO correlations in 7 — 7 production for proton-proton scattering at 
y/s = 5.5 TeV. 



vanish as 4> m - m approaches tt. 

The 7 — 7 correlation functions displayed in Figure 17 indicate that higher order 
corrections do not strongly modify the LO results in this channel either, except near 
the infrared singular point (q T = GeV or z 34 = 1) or in the domain where new 
phase space is available (small m ii or small z 34 ). In the latter case the NLO momen- 
tum configurations modify the LO spectra in a way similar to the 7 — 71° correlations 
(see the spectrum ratios in Figure 18). Although the effect proves tiny, we may also 
remark the lower threshold for the direct process in the transverse momentum dis- 
tribution, now slightly below the 25 GeV cut. However, the most remarkable feature 
when going from LO to NLO in 7 — 7 production deals with the infrared sensitivity 
of observables such as the transverse momentum or the momentum imbalance spec- 
trum. When q T gets small as compared to the diphoton invariant mass, the phase 
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Figure 18: Ratio of NLO correlations over LO correlations for 7 — 7 production. See 
Figure. 17 for details. 

space restriction forces the emitted gluons to be extremely soft. The 5 function 
singularity which appeared in the leading logarithm approximation now spreads in 
phase space, due to the partial cancellation of real and virtual NLO terms, and it is 
broadened at NLO accuracy. Indeed, we notice in Figure 17 the significant correc- 
tions which extend up to roughly q T < 10 GeV. This gives us a typical range in which 
the present perturbative calculation may not be reliable. Since q T /m 11 oc 1 — z w 
such a behavior can also be observed in the momentum imbalance spectrum near the 
singular point, Z34 = 1 where the NLO results start to deviate significantly from the 
LO prediction above z 34 > 0.8. Technically, large terms such as a s In 2 (g 2 /m 2 ) and 
a s In (g^/m^) contribute to the direct process making the present fixed order QCD 
calculation not reliable very near the infrared singular point. For a more accurate 
approximation such large terms should be resummed. Although the one (two) frag- 
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mentation component requires one (two) integration(s) over the scaling variable z 3 
(z 3 and z 4 ), which smears out these large logarithms [7] and make the distributions 
regular at small q T , resummation may affect the shape of the distributions in this 
region. 

Let us now discuss the phenomenological implications of higher order corrections 
to our predictions for the spectrum ratios in lead-lead over those in proton-proton 
collisions. First the moderate higher order corrections (except in specific domains 
for some correlations) clearly indicate that neither the absolute attenuation nor the 
shape of the attenuation in lead- lead collision should be affected too much. In fact, 
the larger presence of gluons at higher order, whose energy loss is stronger than for 
quarks, should be responsible for a slightly more pronounced suppression. In the 
large q T or small z 3A region, where the NLO over LO ratio is the largest due to 
the non-collinear configurations, one can also expect the quenching to prove more 
pronounced. Indeed, such regions were not affected much by the medium in our LO 
prediction as they require the fragmentation of very large k T partons, hence with a 
small energy loss effect. At NLO, however, the large q T and small z 3A domain can 
be reached while keeping the parton energy k T not too large (as compared to u c ), 
provided the relative azimuthal angle between the two particles is small enough. 

We emphasized in the previous section the strong attenuation of diphoton pro- 
duction expected in lead-lead collision near the boundary of phase space, in particular 
at large z 34 and small q T . On the other hand, the presence of the direct process, un- 
affected by parton energy loss, should make the ratio equal to one exactly at z 3i = 1 
and q T = GeV. The competition between the direct and the fragmentation process 
at LO therefore generated discontinuities in the ratio at these specific points (see e.g. 
Figure 13), which should be smoothed at higher order. Based on the present NLO 
calculation in proton-proton reactions, we expect that the quenching of diphoton 
production should start to increase below q T < 10 GeV or z 34 > 0.8. Similarly, the 
discontinuity seen in Figure 13 for 25 GeV photon transverse momenta should be 
smeared as well. 

7. Background 

In the photon transverse momentum range discussed above the background from 
7T° decays will still be appreciable. In order to illustrate this background we briefly 
present here various 7r°— 7r° distributions, using the same asymmetrical cuts as before, 
namely p T ^ > 25 GeV and "p T ^ > 5 GeV (see Figure 19). In this case only the 2f 
mechanism contributes and, as a consequence, the distributions display very similar 
features to the 7 — tt° case. Only the size of the correlations is larger by roughly a 
factor 50. Such distributions should therefore be determined with a great accuracy 
to be subtracted in order to measure the 7 — 7r° and 7 — 7 distributions discussed so 
far. 
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Figure 19: 7r° — 7r° correlations. Both pions are produced at rapidity [-0.5, 0.5] and the 
following cuts are imposed: p T > 25 GeV and p TlT2 > 5 GeV. The conventions for the 
symbols are as in Figure 5. 



Although one may still notice significant effects in nucleus-nucleus collisions, 
medium-modified fragmentation functions cannot be determined through 7r° — 7r° 
correlations in the absence of the If component. Nevertheless, we feel it should be 
interesting to perform such correlations with low p T ^ cuts for both pions which would 
possibly allow to study the spatial distribution of the hot medium. This has been 
achieved for instance by the RHIC experiments who considered the it — tt° azimuthal 
correlations [36]. Theoretically, this would require a complete description of the 
space-time energy density though - available e.g. in hydrodynamical models - which 
go beyond the scope of the present study. Such attempts have been suggested re- 
cently [28, 29]. 
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8. Conclusions 



We have discussed various photon tagged correlations as a tool to study jet frag- 
mentation in the hot medium created in heavy-ion collisions. Correlations functions 
have been computed to leading order in proton-proton collisions at LHC energy. Sim- 
ilar distributions were determined in lead-lead collisions, assuming medium-modified 
fragmentation functions to account for the parton energy loss process in the dense 
medium. 

We show that significant effects could be expected at LHC energy both in the 7 — 
7T° and 7—7 channel. Ideally, the use of asymmetric cuts in the transverse momentum 
of both particles allow the possibility to map out the parton fragmentation functions 
modified by the medium. However several production mechanisms co-exist with 
a relative weight varying with the kinematical variable under consideration. This 
smears somehow the relation between observables and the fragmentation functions. 
Consequently, we found various suppression patterns, depending on the imposed 
kinematical constraints, which should be accessible experimentally. Furthermore, 
the variety of observables presented here should help to constrain the underlying 
model for parton energy loss. 

Although calculations were performed at leading order, NLO corrections have 
also been addressed. In particular, the way higher order corrections could modify 
the expected quenching of 7 — rr° and 7 — 7 spectra is discussed. To be more specific, 
we believe our present LO prediction to be reliable up to roughly z, M ~ 0.8. Finally, 
7r° — 7T° correlation functions were computed so as to give a reference for the expected 
background one could face at the LHC within the kinematic cuts we employed. 
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